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ABSTRACT 

Aims. The Baade-Wesselink method of distance determination is based on the oscillations of pulsating stars. The key 
parameter of this method is the projection factor used to convert the radial velocity into the pulsation velocity. Our 
analysis was aimed at deriving for the first time the projection factor of S Scuti stars, using high-resolution spectra of 
the high-amplitude pulsator AI Vel and of the fast rotator f3 Cas. 

Methods. The geometric component of the projection factor (i.e. po) was calculated using a limb-darkening model of 
the intensity distribution for AI Vel, and a fast-rotator model for f3 Cas. Then, using SOPHIE/OHP data for /3 Cas and 
HARPS/ESO data for AI Vel, we compared the radial velocity curves of several spectral lines forming at different levels 
in the atmosphere and derived the velocity gradient associated to the spectral-line-forming regions in the atmosphere 
of the star. This velocity gradient was used to derive a dynamical projection factor p. 

Results. We find a flat velocity gradient for both stars and finally p = po = 1-44 for AI Vel and p = po = 1.41 for /3Cas. 
By comparing Cepheids and 8 Scuti stars, these results bring valuable insights into the dynamical structure of pulsating 
star atmospheres. They suggest that the period-projection factor relation derived for Cepheids is also applicable to 

5 Scuti stars pulsating in a dominant radial mode. 
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1. Introduction 

Determining distances in the Universe is not a trivial task. 
From our Galaxy to the Virgo Cluster, distances can be de- 
rived using the peri od-luminosity relation (PL) of Cepheids 
(jRiess et al.ll2009aH bT). However, this relation has to be cal- 
ibrated, using the Baade- Wesselink method o f distance de- 
termination for instance (jStorm et al.l 120 lTal fbl) . The prin- 
ciple of this method is simple: after determining the angu- 
lar diameter and the linear radius variations of the star, 
the distance is derived by a simple ratio. Angular di- 
ameter variations ca n be measured using interferometry 
dKervella et al. 20041) or the in frared surface brightness re- 
lation ( Gieren et al ri998ll2005lh The linear radius variation 
is measured by integrating the pulsation velocity (hereafter 
^puis) over one pulsating cycle. However, from observations 
we have only access to the radial velocity (VJ- a d) because of 
the projection along the line-of-sight. The projection fac- 
tor, used to convert the radial velocity into the pulsation 
velocity, is defined by p = V vu \ a /V ra ^. There arc in prin- 
ciple three sub-concepts involved in the Baade-Wesselink 



projection factor: (1) the geometric projection factor po, 
which is directly related to the limb-darkening of the star 
(see Sect. [3]), (2) the correction / gra d due to the velocity gra- 
dient between the spectral-line-forming region and the pho- 
tosphere of the star; this quantity can be derived directly 
from observations by comparing different lines forming at 
different levels in the atmosphere (see Sect. 2 and|4|), and 
(3) the correction / - g due to the relative motion between 
the optical and gas layers associated to the photosphere 
(see Sect. 2]). For a detailed analysis of the p- factor decom- 
position we refer to iNardetto et al.l (|2007f ) . The projection 
factor is then defined by p — Po/grad/o-g- In the follow- 
ing, we apply this decomposition of the projection factor 
(originally developed for Cepheids) to the 8 Scuti stars AI 
Vel and j3 Cas. The impact of non-radial modes of 6 Scuti 
stars on the project ion factor is a very difficult question 
studied previously (iDziembowski _ 1977 j Balona fc Stobid 



* This work uses observations made with the HARPS instru- 
ment at the 3.6 m telescope (La Silla, Chile) in the framework 
of the LP185.D-0056 and with the SOPHIE instrument at OHP 
(France) . 



Il97f* IStamford fc Watsonlll981t lHatzesHl99fih . The nartic 
ular cases of AI Vel and /3Cas are discussed in the con- 
clusion. This paper is part of the international Araucaria 
Project, whose purpose is to provide an improved local cali- 
bration of the extrag alactic distance scal es out to distances 
of a few megaparsecs (|Gieren et al.ll200"5l ). In this context, S 
Scuti stars are extremely interesting since it has been shown 
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recently that they follow a PL relation (jMcNamara et all 
120071: IPoretti et aUl2008lh 

2. Spectroscopic observations of 5 Scuti stars 

AI Vcl (HD 69213, A9IV/V) is one of the most often 
observed double-mode, high-amplitude 8 Scuti stars. This 
star pulsates in the fundamental and first overtone radial 
modes with a well-c onstrained period ratio Pi/Po of 0.77 
(jPoretti et ai ] |2005D. In addition t o P n = 0.111574d and 
Pi = 0.0862073 d. IWalraven et al.l (|1992T) clearly detected 
two other periods, tentatively identified as the third and 
fifth radial overtones. We observed AI Vel using the HARPS 
spectrograph mounted at the ESO 3.6-m telescope. We 
analysed 26 high signal-to-noise ratio (S/N~140) spectra 
taken in the high-efficiency mode (EGGS, R=80 000) on 
the night of January 9-10, 2011. We identified 53 metal- 
lic unblended spectral lines (ranging from 3 780 to 6 910 A) 
relevant for the determination of radial velocities. Figure Q] 
{left panel) shows the behaviour of the mean spectral line 
profile along the pulsation phase. The shifts due to the ra- 
dial modes clearly dominate the line profile variations. 
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Fig. 1. Mean profiles of the HARPS spectra of AI Vel (left 
panel, T =JD 2443176.00) and of the SOPHIE spectra 
of /3 Cas (right panel, T =JD 2438911.88). The pulsation 
phase is given on the y-axis. A strong broadening is clearly 
seen for /3 Cas due to its high rotation velocity. We com- 
puted the mean line profiles of AI Vel and (3 Cas spectra by 
means of a deconvo lution process using the LSD software 
(|Donati et al.lll997l) . 

The target (3 Cas (HD432, F2III/IV) is a low- 
amplitude S Scuti star. We observed (3 Cas with the 
SOPHIE instrument (R = 75 000) at the OHP 1.93-m tele- 
scope on the night of September 30, 2011. We collected 241 
high- resolution spectra with a mean S/N of 100. We could 
distinguish only height unblended spectral lines relevant for 
the spectral analysis beca use of the strong rota tional broad- 
ening (Fig. [I] right vanel) . iRiboni et al.l (|l994l ) showed that 



the star is a monoperiodic pulsator at the detection limit of 
ground-based photometric measurements, with a pulsation 
period o f P = 0.101036676 d. The mode identifi cation is 
unclear (|Rodriguez et al.lll992l Ikiborii et al.lll994h . Today, 
the distance of (3 Cas is known to b e 16.8 pc from the 
Hipparcos parallax ([van Leeuwenll2007l) . Thus we could ob- 
tain the absolute magnitude Af^=1.14 f rom the apparent 
magnitude V=2.27. The PL relations (jMcNamara et all 
l2007t IPoretti et al.ll2008l) supply a fundamental radial pe- 
riod of about 0.15 d at this My value. Therefore, the ob- 
served period is similar to that expected for the second 
radial overtone. We attempted a mode identification from 
our spectroscopic data using the FAMIAS3 software, bmce 
we dealt wit h a fast rot ator, we used the Fourier parameter 
fit method (|Zimall2006lh Imposing the frequency 1/P, the 
results from spectroscopy point towards an axisymmetric 
mode, without a clear indication on the I- value. Since (3 Cas 
is seen almost pole-on (i = 19.9 ± 1.9°. IChe et~aTI (|2011l )). a 
low-degree, axisymmetric, nonradial mode mimics the pul- 
sation behaviour of a radial mode very well. On the basis of 
these considerations, we treated j3 Cas as a radial monope- 
riodic pulsator for our purposes. We also used the mean 
line profiles of f3 Cas to estimate the v sin i values f rom the 
posit ion of the first zero of their Fourier transforms (|Carrolll 
Il933f ). This approach is possible only for objects where the 
rotational broadening is dominant with respect to the other 
sources of line broadening (e.g., instrumental effect, micro- 
turbulence), which is always the case for (3 Cas, but not 
for AI Vel, where we were unable to use the Fourier trans- 
form method on the narrower lines (vsini < 10 km s _1 ). 
The radial velocity values of the observed profiles of (3 Cas 
range from 5.3 to 11.6 km s _1 and the vsini values from 
74.0 to 77.5 km s _1 (Fig[TJ right panel). We could also 
determine a mean v sin i value from the average of all the 
mean profiles and obtained 75.72±0 .14 km s -1 . This value 
is consistent with literature values dBernacca fc Perinottol 
119701: lUesugi fc Fukudalll970t ISchroder et al.ll2009h . When 
considering i = 19.9°, this means that /3Cas is an intrinsic 
fast rotator, wi th a velociy of v ro t ~220 km s _1 which is 
consistent with lChe et al.l (|201lD . 

Finally, for both stars, the centroid radial velocity RV C 
(or the first-moment radial velocity) and the line depth D 
are derived as a function of the pulsation phase for each 
selected spectral line. These data are used in Sect. 4. 

3. The geometric projection factor po 

Considering a limb-darkened pulsating star in rotation with 
a one-layer atmosphere, the projection factor is purely ge- 
ometric. Thus, p = = po. The radial velocity is then 
defined by 
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where Dr is the surface of the stellar disc of radius R, 
and I(x, y, A) the limb-darkened continuum intensity dis- 
tribution considered at the wavelength of observation A 
defined by I(x,y,X) = I (l - u x + u^l - (x 2 + y 2 )), 
where u\ is the linear limb-darkening coefficient from 



1 Developed in the frame work of the FP6 European 
Coordination Action HELAS (http://www.helas-eu.org/!. 



2 



G. Guiglion and collaborators: The Baade-Wesselink projection factor of the 8 Scuti stars 



IClaret fc Bloemenl (|201lh . Considering T cff = 7400K and 
log g = 3.5, we find Mr = 474 ± 0.025 in the i?-band 
from IClaret fc Bloemenl (|2011[) . Using Eq. 1, we deduce 
a value of the geometric projection factor for AI Vel of 
Po = 1.43 ± 0.0 1. jPp is assumed to b e constant with the 
pulsation phase (jNardetto et alj|2~004l ). This value is higher 
than what we generally obtain for Cepheids (typically 1.37 
to 1.41, see Figf2](fop panel)). 

The geometrical shape of ft Cas and its intensity distri- 
bution arc distorted by its high rotation rate. The geomet- 
ric projection factor depends on the inclination of the star's 
rotation axis compared to the line-of-sight. If the star's ro- 
tation axis is along the line-of-sight (i = 0°), the star is 
observed pole-on and is seen as a disc. For i > 0° the star 
has an ellipsoidal shape. 
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where a and b are the semi-major and semi-minor axis of 
the ellipse. In Fig. [2] (bottom panel), we show the geomet- 
ric projection factor (po) as a function of i. Fig. [3] presents 
the modelled intensity distributions for several inclinations. 
This relation is extremely interesting because it shows that 
the inclination of a fast-rotating star can have an impact 
of more than 10% on the projection factor. Of course, it 
also depends on the rotation velocity of the star: the higher 
the rotation velocity (for a given inclination), th e lower the 
projec tion factor. Using the inclination found bv lChe et al.l 
(|201ll ). i = 19.9 ± 1.9°, we finally find a geometric projec- 
tion factor for (3 Cas of po = 1.41 ± 0.02 (averaged over the 
three wavelengths considered). 

Fig. 3. Modelled intensity distributions of (3 Cas at 
6 000 A for i = 5° (left), i = 50° (middle), i = 90° (right). 
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Fig. 2. Top: po as a function of the limb-darkening param- 
eter u\. The red box indicates the uncertainty on po for 
the S Scuti AI Vel. The blue box indicates the typical val- 
ues of u\ and po for C epheids. The dots co rresponds to 
the relation provided bv lNardetto et all (|2006|) . Bottom: po 
as a function of the inclination of the fast rotating star 
/3Cas for three different wavelengths (A = 6 000 A (□), 
A = 6500 A (A), and A = 7000 A (x)). The red box in- 
dicates the uncertainty on po for (3 Cas. The case of an 
uniform clonged disc is over-plotted (+), and we find that 
Po = 1.5 for i — 0°, as expected for a circular uniform disc. 



Using the fundament al paramet e rs of the modified von 
Ziepel model fou nd by Cheetah (12011) and the rotat - 
ing stars model bv lDomiciano de Souza et all (|2002l . 
we derived the intensity distribution in the continuum 
for different inclinations of the star (from i = 0° to 
i = 90° with a step of 5°) and for three wavelengths: 
A = 6 000, 6 500 and 7 000 A. Using these intensity maps, we 
can easily calculate the geometric projection factor. Indeed, 
for an ellipsoid, U ra d is then defined by 



4. Dynamical structure of the atmosphere 

By comparing the 2K-amplitude (defined as the ampli- 
tude of the RV C curve, hereafter Ai?U c ) with the depth 
(D) of the 53 spectral lines selected in the case of AI Vel, 
one can directly measure the atmospheric velocity gra- 
dient in the part of the atmosphere where the spectral 
lines are formed (Figure El top- left). To quantify the im- 
pact of velocity gradient on the projection (/ gra d), we do 
not need to derive the velocity gradient over the whole 
atmophere, but only at the location of the forming re- 
gions of the spectral lines used to derive the distance of 
the star. We therefore performed a linear regression ac- 
cording to the relation ARV C = ao D + bo- We obtained 
Ai?U c = [-0.40 ± 0.53]D + [32.87 ± 0.23] kms" 1 (Fig. H 
top- right). I n principle, / gra d dep ends on the spectral line 
considered (jNardetto et all 120071 ): / grad = bo/(a D + b ). 
Here, we find that / gra d is typically the same for all spec- 
tral lines (/grad = 1.01 ± 0.01), which is consistent with no 
correction of the projection factor due to the velocity gra- 
dient. The uncertainty on / gra d is derived from the errors 
on ao and bo- 

Figure|4] (bottom left) presents the interpoled RV C curve 
of ,3 Cas for the Fel spectral line (A = 4 508.288 A). We 
clearly see an increase of the amplitude of the radial ve- 
locity curve (~ 4.6 ± 0.9% per cycle). Moreover, the ra- 
dial velocity curves have several minima and maxima and 
we can easily deduce a period of pulsation. We find P = 
0.10046 ± 0054 d . Our value agrees well with that of 
iRiboni et all (| 19941) . 

In Fig. 0] (bottom right), ARV C is plotted as a function of D 
for cycle 2 (see the vertical line in the figure). For /3 Cas, the 
range of the spectral line depth is eight times lower com- 
pared to AI Vel. The velocity gradient is / gra d = 0.64±0.82. 
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Fig. 4. Top: RV C as a function of the MJD in the case 
of the Fell 5 234.625 A spectral line, and amplitude of the 
RV C curves as a function of the spectral line depth (D) 
for AI Vel. Typical error bars are indicated in the upper 
right of each panel. Bottom: The same for (3 Cas f the Fell 
4 508.288 A spectral line. 



The large error bar prevents any physical discussion about 
the exact value of / gra d- We note that / gra d =0.64 implies 
a huge velocity gradient in the star's atmosphere, which 
seems unrealistic. Since it is consistent with our determi- 
nation, we assumed / gra d = 1-0. Considering ciq = 0, we 
obtained 6 = 4.07 ± 0.25 with a reduced \ 2 01 1-0, which 
provides an uncertainty on / gra d of 0.14. The / _ g correc- 
tion, which is the last component of the projection factor 
decomposition, cannot be measured from observations. To 
estimate the differential velocity between the optical and 
gas layers at the photosphere of the star, we need a hydro- 
dynamic model. 

However, modelling the pulsating atmosphere of 5 Scuti 
stars is not an easy task because of (1) cyclc-to-cycle vari- 
ations (non-radial modes) and (2) fast rotation in some 
cases (as for /3Cas). Ho wever, / n _„ have been s tudied inten - 
sively in the Cepheids (INardetto et alj|2004 120071 120111) . 
and it seems that there is a linear relation between / Q _ g 
and logP: / „ g = [-0.023 ± 0.005] logP + [0.979 ± 0.005]. 
Moreover, we have a theoretical value of / _ g for the short- 
period /3-Cephei aLup (P = 0.2598 d), / Q _ g = 0.99 ± 0.01, 
which seems to b e consistent with th e logP-/ - g rela- 
tion of Cepheids (|Nardetto et al.l I2012D . For our study, 
we propose to extend this law for the 5 Scuti (3 Cas and 
AI Vel. Considering P = 0.11157 d for AI Vel and P = 
0.10046 ± 0.00054 d for /3Cas (from this paper), we find 
/ _ g = 1.00 ± 0.02 (0.02 is a conservative arbitrary un- 
certainty) for both stars, which basically means no pho- 
tospheric correction to the projection factors. 



5. Discussion 

We can now calculate the projection factor p, using the 
relation p — _po/ g r a d/o- g - We find p = 1.44 ± 0.05 for AI 



Vel and p = 1.41 ± 0.25 for /3 Cas. However, the generali- 
sation of this study to any 5 Scuti stars is presently lim- 
ited since we have to study the impact of multi-modes, 
in particular non-radial ones, on the projection factor. 
This complicated effect has been studied by several au- 
thors for the bisector method o f the radial velocity de- 



termination (iDziembowskl Il977t iBalona k Stobiel f 19791: 



iStamford k Wat son! 1 19811 : lHatzeE 1 1996T ) . We assume a star 



pulsating in two modes, one radial and one non-radial. The 
projection factor can be defined as p — a x p T + a nr Pnr, with 
a r and a m the relative contributions of the velocity am- 
plitudes of the radial and non-radial modes to the pulsa- 
tion (a r + a m = 1). p T is our previous decomposition of 
the projection factor in the case of a purely radial mode 
Pr = Po/gr a d/o- g , while p nl is the projection factor in the 
case of a purely non-radial mode. Using Eq. 4 of lHatzesI 
(|l996l ). we find that p nr = poe +ke , where k = 0.15 in the 
case of a non radial p-mode and k = 1.2 in the case of 
a non-radial g-mode. I is the spherical harmonic degree. 
This relation is derived for the first moment (i.e. the ra- 
dial velocity) determination only, which is independent of 
the star's rotation. This equation can be used when £ = m 
only (where m is the spherical harmonic order). Additional 
work is necessary to derive it when I ^ m. We emphasize 
also that the higher £(= m), the higher is the non-radial 
projection factor. This is expected since for high values of 
t (and this is qualitatively also true when £ ^ m) there are 
more red- and blue-shifted velocity zones on the star that 
cancel each other in the integrated line profile, which leads 
to a lower amplitude of the non-radial velocity curve, and 
in turn a hi gh valu e of th e non-radial projection factor (see 
Eq. 1&2 of lHatzesI (I1996D V 

Wc performed our study on AI Vel and (3 Cas under the 
assumption of monoperiodic radial pulsation. Our results 
(without non-radial correction) are consistent (at the la 
level) with the period-projection-factor (P p) relation p = 
-0.07 1 ± 0.020] logP + [1.311 ± 0.019] bv lLanev k Jonen 
2009) applied for classical and dwarf Cepheids (it corre- 
sponds to p = 1.38 ± 0.02 for AI Vel and j3 Cas). To derive 
these values, Laney k Joner simply compared the distance 
of the stars obtained from the PL relation with the dis- 
tances from the photometric version of the Baade-Wesselink 
method. This suggests that the eventual non-radial compo- 
nents of f3 Cas have probably a negligeable impact on the 
projection factors (which means a nr ~ 0). Interestingly, 
if we use the Pp relation obtained for classical Cepheids 
by iNardetto et all (120071) . p = [-0.064 ± 0.020] logP + 
[1.376 ± 0.023] (derived with the first moment mcthocQ) 
to derive the projection factors of the two 8 Scuti stars, 
we find p = 1.44 ± 0.01, which is consistent with our val- 
ues. This seem s to sh ow that the Pp relation provided by 
INardetto et alj (|2007l ) for single lines is also applicable to 
5 Scuti stars pulsating in a dominant radial mode. In addi- 
tion, and as already shown, for fast- rotating S Scuti stars, 
an intrinsic dispersion of the Pp relation due to the random 
orientation of the rotation axis has to be considered. 
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2 We remind that using t he cross-correlat i on me thod, one has 
to use the Pp relation from INardetto et al.l l|2009f ). 
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